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ABSTRACT 

We measure the average mass properties of a sample of 41 strong gravitational lenses at moderate 
rcdshift (z ~ 0.4 - 0.9), and present the lens redshift for 6 of these galaxies for the first time. Using the 
techniques of strong and weak gravitational lensing on archival data obtained from the Hubble Space 
Telescope, we determine that the average mass overdensity profile of the lenses can be fit with a power- 
law profile (AS oc ^-°- 86 ±o.i6-j -j- na ^ jg w jthin l-er of an isothermal profile (AE cx i?" 1 ) with velocity 
dispersion a v = 260 ± 20 km s _1 . Additionally, we use a two-component dc Vaucouleurs+NFW 
model to disentangle the total mass profile into separate luminous and dark matter components, and 
determine the relative fraction of each component. We measure the average rest frame V-band stellar 
mass-to-light ratio (Ty = 4.0±0.6 h Mq/Lq) and virial mass-to-light ratio (ry = 300±90 h Mq/Lq) 
for our sample, resulting in a virial-to-stellar mass ratio of M vir /M^ = 75 ± 25. Relaxing the NFW 
assumption, we estimate that changing the inner slope of the dark matter profile by ~20% yields 
a ^30% change in stellar mass-to-light ratio. Finally, we compare our results to a previous study 
using low redshift lenses, to understand how galaxy mass profiles evolve over time. We investigate the 
evolution of M vir /M ir {z) = a(l + z)' 3 , and find best fit parameters of a = 51 ± 36 and /? = 0.9 ± 1.8, 
constraining the growth of virial to stellar mass ratio over the last ^7 Gigayears. We note that, by 
using a sample of strong lenses, we are able to constrain the growth of M V i r /M*(z) without making 
any assumptions about the IMF of the stellar population. 

Subject headings: dark matter - gravitational lensing - galaxies: elliptical and lenticular, cD - galaxies: 
evolution - galaxies: structure 



1. INTRODUCTION 

Observations over the last few decades have sug- 
gested that galaxies are embed ded in an extended , 
diffuse halo of dar k matter (e.g. iPeterson et al.1 119781: 
White fc Reeslll978l : [Rubin et al.lll979tlBlumenthal et al l 
19841) . whic h contributes up to 95 percent of the 
total mass dHoekstra et al.1 [20051 Uiang fc Kochanekl 
20071: lAvila-Reese et all 120081 ). Because of the over- 
whelming fraction of dark matter, numerical simu- 
lations dealing with the formation and assembly of 
galaxy mass often focus on dark matter alone (e.g., 
the Via Lactea Simulation; IKuhlen et al.1 I2008D . and 
these simulations have made very specific predictions 
about the overall profile of galaxy-sized mass distribu- 
tions (e.g.. iNavarro et all 119971: iMoore et all 119981: Uing 
20001: IStoehrl 120061: IDiemand et al .1 120071 : ISchmidt et al 
20081 ). However, when compared to observations, these 
simulations often fail to accurately rec reate the ob- 
served shape of galaxy mass profiles (|Saluccil 12001 
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Gentile et all 120041: ISimon et all [200l IGentile et all 
20071: IKuzio de Narav et all l2008MSeigar et all 120081: 
Sand et al. I 1200811 This disagreement is attributed 



largely to the complex physics of baryon interaction (e.g. 
frictional dissipation and radiative cooling) which is still 
not well understood and is difficult to model. 

There are many possible methods capable of measur- 
ing the distribution of matter on small (< 1 Mpc) scales, 
and techniques used for characterizing a radial mass 
profile include measuring the rotation curves from plane- 
tary ne bulae (iRomanowskv et alll2003l: lArnaboldi et all 



2001 IMerrett et all 120061: Ide Lorenzi et all 120081: 
Napolitano et al. 20091 ), kin e matics of stellar pop- 



ulations dBertin et all 119941 : iCappellari et all 120061 : 
Ivan der Wei fc van der Marell 120081) and H I gas 
(lUson fc Matthews! 120031: I Jackson et all 120041: 
lAndersen et all 120061: iMatthews fc Usonl 1200811 or 



the temperature of X-ray e mitting gas ( Humphrey et all 
l2006l : IChurazov et al.ll20M) . 

In this paper, we use a powerful al ternative tech- 
nique : gravitationa l lensin g (e.g.. | Brainerd et all 
19961: iFischer et all [20001: IWilson et all [200iT 
Kleinheinrich et all 120061: iParker et all 120071) . The 



power of lensing is due to the fact that the technique 
is able to directly trace the total mass (luminous + 
non-luminous) enclosed within a given radius without 
needing to make any assumptions about the dynamical 
state of the mass in question. In addition, lensing does 
not rely on the presence of kinematic tracers, which are 
often only visible in very low-redshift galaxies, and even 
then are typically not present in the outer halo regions 
where dark matter dominates. 
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Furthermore, since we analyze a sample of strong 
lenses, we are able to combine the mass constraints deter- 
mined from both strong lensing and weak lensing. This 
allows us to probe the total mass distribution over a 
much wider range of physical distances than by using 
weak lensing alone, and also allows us to constrain the 
properties of stellar mass distribution without selecting 
a specific stellar initial mass function (IMF). Instead, we 
rely on directly-derived properties: lensing-inferred mass 
and total luminosity (although we do make an assump- 
tion about the shape of the dark matter halo). 

For early-type galaxies, several lensing-based studies 
have shown the typical shape of the mass density pro- 
file to be consistent with an isothermal (p(r) oc r~ 2 ) 
model, given the uncertainties on the mea surements, be- 
tween distances of ~ 50 — 30 ft -1 kpc (jSheldon et al.l 
12004 IMandelbaum et"aTll2006D . However, due to the low 
signal-to-noise ratio (SNR) on the lensing signal obtained 
from an individual galaxy, these studies rely on stacking 
a large number of galaxies into a single sample to de- 
termine the average profile of the entire stack. For this 
work, we utilize deep imaging from the Hubble Space 
Telescope (HST) . In general, space- based data will have 
a much higher density of background galaxies (n s ) when 
compared to data taken from the ground over an iden- 
tical exposure time. This increased background density 
count not only enables us to detect a weak lensing sig- 
nal much closer to the center of the lensing galaxy, but 
also allows for a significant detection with fewer stacked 
galaxies (Ni ens ), because the expected weak lensing SNR 
scales as VNicns n s . The ability to work on a small sam- 
ple size is particularly useful, given the relative paucity 
of known strong lenses at all redshifts. 

In this wor k , we are able to extend the results of 
iGavazzi et al.1 (j2007l hereafter G07), who conducted a 
similar mass profile study on a small (22 lens) sam- 
ple of low redshift (z ~ 0.2) early- type strong gravita- 
tional lenses colle cted from the Sloa n Lens ACS Survey 
(SLACS; see e.g. iBolton et al1[200l [20081 ) over a wide 
range of physical distances (~ 3 — 300 h~ x kpc). From 
this work, it was shown that, like their non lensing coun- 
terparts, the total mass density profile of strong lens- 
ing ellipticals could be described by a roughly isothermal 
model (although they note that the comparison between 
the data and the model was done without a formal fit). 
This was thought to be due to the combination of bary- 
onic and dark matter, the so-called Bulge-Halo Conspir- 
acy, in which a relatively steep luminous matter profile 
and a shallower dark matter profile combine in such a way 
that the total mass distribution of the galaxy can be we ll 
approximated by an isothermal model (|Treu et"aT1l2006t ). 
By comparing this data set to our moderate redshift sam- 
ple, we thus explore the evolutionary trends of massive 
early-type galaxies. 

This paper is organized as follows: In Section 2, we de- 
scribe our lens sample, the lens selection criteria, and the 
data reduction pipeline used for this work. At this point, 
the interested reader may turn to the appendix, where we 
describe our techniques for analyzing weak lensing shear, 
paying specific attention to signal to noise optimization 
and the removal of systematic errors. In Section 3, we 
use a two-component model to determine the contribu- 
tions of the stellar and dark matter components to the 



total profile. We present all of our results in Section 4, 
the reader interested in the science may skip directly to 
this section. We discuss the results and compare them to 
the SLACS results in Section 5, looking for evolutionary 
differences between the samples. Finally, we summarize 
our results in Section 6. 

Throughout this work, we assume a flat cosmology 
with H = 100 h km s" 1 Mpc" 1 , Ct M = 0.3, and 
Q\ = 0.7. All magnitudes presented in this paper are 
AB magnitudes. 

2. DATA 

2.1. Imaging Data 

For this analysis, we focus on a sample of 41 strong 
gravitational lenses at moderate redshift (z ~ 0.4 — 
0.9, ^median ~ 0.6). A majority of the sample was 
observed as part of the CASTLES prograrrQ, but 
it als o contains lenses found in the COSMOS sur- 
ye^ dFaure et alj|2008t) and the Extended Groth Strip 
(jMoustakas et al.ll2007t) . as well as targeted exposures 
of individual lens systems. Unlike other data sets that 
have been selected in some uniform way (such as the 
SLACS lenses), our lens sample is drawn from a variety 
of sources and selected simply for being lenses with early- 
type morphologies in the desired redshift range. How- 
ever, we do note that systems known to be strongly af- 
fected by the presence of a galaxy cluster (Q0957+561, 
SDSS1004+4112, and B2108+213) were excluded from 
the sample. A full list of the lenses can be seen in Ta- 
bled! 

Each system in the sample was chosen based on the 
existence of publicly available HST images, subject to 
the following criteria: the system needed to be observed 
using the Wide Field Chan nel (WFC) of the Advanced 
Camera for Surveys (ACS; iFord et all 119981 ). and the 
total exposure needed to be at least ~ 1800 seconds 
through the F814W filter. This exposure time yields 
background densities of ~ 70 galaxies arcmin -2 or more. 

The imaging data for all the lenses were obtained 
through the Multimission Archive at STScI (MAST0), 
from programs G0-9744 (CASTLES: PI Kochanek), GO- 
9822 and 10092 (COSMOS: PI Scoville), GO-10134 
(EGS: PI Davis), GO-10158 (B1608+656: PI Fass- 
nacht), and GO-9450 (B0218+357: PI Jackson). The 
lens systems that had been specifically targeted by HST 
(B1608+656, B0218+357, and the CASTLES lenses) are 
located at the WF1 target point, whereas the serendip- 
itously discovered COSMOS and EGS lenses appear at 
random positions in the ACS field of view. 

In addition to the F814W data, some lens systems have 
secondary imaging in either the F555W or F606W filter. 
These extra data could prove to be useful in future re- 
finements to this study, especially in using color selection 
to reject foreground interlopers from the population of 
background sources (see the appendix). However, we do 
not use the information in our present weak lensing anal- 
ysis: a significant fraction of our sample has only F814W 
imaging, and for systems that do have multi-band data, 
the redshift /color relationship can be highly degenerate 

7 http:/ /cfa-www. harvard.edu/castles/ 

8 http:/ /cosmos. astro. caltech.edu/ 

9 http://archive.stsci.edu/ 
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Table 1 

Lens System Data. 



Lens Name Program Exp. Time Z; a,b z s h D\ iy(z;) c S cr ; t References 





ID 


(s) 






(h^ 1 Mpc) 




(M pc" 2 ) 










B0218+357 


9450 


4320.0 


0.68 


0.96 


1020 


0.320 


5097 [11, [21, 


3 


SDSS0246-0825 


9744 


2288.0 


0.724 


1.68 


1045 


0.293 


5035 








1 


CFRS03P1077 


9744 


2296.0 


0.938 


2.941 


1137 


0.188 


8965 








5 


HE0435-1223 


9744 


1445.0 


0.454 


1.689 


842 


0.486 


3702 


[61. 


7 


B0445+128 


9744 


5228.0 


0.557 




931 


0.406 


3969 








8 


B0631+519 


9744 


2446.0 


0.62 




979 


0.360 


4267 








9 


J0816+5003 


9744 


2440.0 






825 


0.500 


3678 






10 


B0850+054 


9744 


2296.0 


0.59 


3.93 


957 


0.381 


4112 






11 


SDSS0903+5028 


9744 


2444.0 


0.388 


3.605 


761 


0.552 


3645 






12 


SDSS0924+0219 


9744 


2296.0 


0.394 


1.524 


763 


0.550 


3644 


[131, 


11 


J1004+1229 


9744 


2296.0 


0.95 


2.65 


1141 


0.183 


9397 






15 


HE1113-0641 


9744 


1062.0 


0.75 


1.235 


1059 


0.278 


5298 






16 


Q1131-1231 


9744 


1980.0 


0.295 


0.658 


635 


0.646 


3803 






17 


SDSS1138+0314 


9744 


2296.0 


0.45 


2.44 


831 


0.495 


3685 


[181, 


19 


SDSS1155+6346 


9744 


1748.0 


0.176 


2.89 


430 


0.780 


4782 






20 


SDSS1226-0006 


9744 


2296.0 


0.52 


1.12 


899 


0.435 


3841 


[18], 


19 


B1608+656 


10158 


9744.0 


0.63 


1.39 


987 


0.353 


4774 


21], 


22 


WFI2033-4723 


9744 


2085.0 


0.66 


1.66 


1007 


0.333 


4516 


[18], 


23 


COSMOS5857+5949 


9822 


2028.0 


0.39 




763 


0.550 


3646 






21 


COSMOS5914+1219 


9822 


2028.0 


1.13 




1169 


0.148 


15755 






21 


COSMOS5921+0638 


9822 


2028.0 


0.551 


3.15 


926 


0.411 


3948 


[24], 


25 


COSMOS5941+3628 


9822 


2028.0 


0.88 




1124 


0.204 


7834 






24 


COSMOS5947+4752 


10092 


2028.0 


0.345 




706 


0.595 


3681 


24 




26 


COSMOS0012+2015 


9822 


2028.0 


0.378 




749 


0.562 


3649 


21 




26 


COSMOS0013+2249 


9822 


2028.0 


0.346 




707 


0.594 


3680 


24 




26 


COSMOS0018+3845 


9822 


2028.0 


0.71 




1037 


0.301 


4908 






21 


COSMOS0038+4133 


10092 


2028.0 


0.738 




1121 


0.208 


7583 [241, [29 


COSMOS0047+5023 


9822 


2028.0 


0.87 




1105 


0.226 


6729 






21 


COSMOS0049+5128 


9822 


2028.0 


0.337 




695 


0.603 


3694 


"24], [26 


COSMOS0050+4901 


10092 


2028.0 


0.960 




1144 


0.179 


9795 


"241, 26 


COSMOS0056+1226 


9822 


2028.0 


0.361 


0.81 


727 


0.579 


3661 [241, [261, [29 


COSMOS0124+5121 


9822 


2028.0 


0.84 




1101 


0.231 


6546 






21 


COSMOS0211+1139 


9822 


2028.0 


0.920 




1124 


0.201 


7834 [24], [29 


COSMOS0216+2955 


9822 


2028.0 


0.608 




1013 


0.326 


4588 24,29 


COSMOS0227+0451 


10092 


2028.0 


0.89 




1121 


0.208 


7583 






21 


COSMOS0254+1430 


10092 


2028.0 


0.417 


0.779 


825 


0.500 


3678 


24 


) 


29 


J095930.93+023427.7 


9822 


2028.0 


0.892 




1122 


0.207 


7150 


27 


! 


29 


J100140.12+020040.9 


9822 


2028.0 


0.879 




1117 


0.213 


6987 


27 


■ 


29 


"Anchor" 


10134 


2100.0 


0.463 




845 


0.483 


3707 






28 


"Cross" 


10134 


2100.0 


0.810 


3.40 


1088 


0.246 


6060 






28 


"Dewdrop" 


10134 


2100.0 


0.580 


0.982 


950 


0.389 


4068 






28 



a Lenses with no known z% were placed at a fiducial rcdshift z — 0.6. 

k For the COSMOS lenses, redshift values with 3 significant figures were determined spectroscopically, and redshift values with 2 significant figures 

are determined photometrically, according to the catalog of Ilbcrt ct al. (2009). 

c w(zi ) is the weight value associated with each lens field, described in the appendix. 

R eferences: [1]: IPatnaik et al.l 119931). 12]: IBrowne et al.l |Tl993l). [3]: [Lawrence! fl996D. 141 : llnada et al.l 120051). [5] : ICrampton et al.l (I2002f). 
[6]: IWisotzki et al. (2002), [7]: Morgan et al. (2005), [8 ]: lArgo eT al. (2003), [9]: York ct al. (2005), [10]: Lchar et al. (2001), [11]: Big gs et all 
(2003), [12]: Johnston ct al. (2003), [13]: Inada ct al. (2003), [14]: Eigenbrod ct al. (2006. a), [15]: Lacy ct al. (2002), [16]: Blackburnc ct al. 
120081), [17]: Slusc ct al. (20®, [ 18]: [E igenbrod ct al. (2006, b), [19]: Inada ct al. (2008). [20]: Pindor ct al. (2004j), [21]: IMvers et aTl U995I) . 
[22]: [jassiiacht ct al. ( 1996), [23]: Morgan ct al. (2004), [24]: Faurc ct al. (200|), [25]: [Anguita ct al. (200[|), [26]: C. Faure et al. in preparation, 
[27]: Jackson (200Sl), [28] : Moustakas ct al. (2007]), [29]: This paper 



when dealing with only two filters. Thus, we do not ex- 
pect significant improvement in the signal. 

The raw data were processed through a reduction 
pipeline created for the HST Archive Galaxy-scale Grav- 
itational Lens Survey (HAGGLeS, P.J. Marshall et al., 
in preperation) , which we briefly describe. First, each 
individual raw exposure of a given lens system is cal- 
ibrated, using the calacs package in STSDASFI. a soft- 
ware system built on top of IRAF0. Once calibrated, the 

10 STSDAS is a product of the Space Telescope Science Institute, 
which is operated by AURA for NASA 

11 IRAF (Image Reduction and Analysis Facility) is distributed 
by the National Optical Astronomy Observatory, which is operated 
by the Association of Universities for Research in Astronomy, Inc., 
under cooperative agreement with the National Science Founda- 



exposures are aligned and combined into a single stacked 
image using multidrizzlc (Ko ckemoer et aLll2002[ ). Since 
even slight misalignments can introduce systematic er- 
rors in the weak lensing signal, the HAGGLeS pipeline 
refines the image alignment derived from the astrometric 
header by cross-correlating the positions of bright, well- 
defined objects (bright unsaturated stars, nebular knots 
in spiral galaxies, etc.) in each exposure in order to look 
for residual shift or rotation misalignments. These resid- 
ual shifts are then fed to multidrizzlc, and this process 
is repeated until the shift refinements become negligibly 
small. After final alignment and combination, the com- 
posite output image is resampled in multidrizzlc from the 

tion. 
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Table 2 

Relevant SExtractor Parameters. 



z\ y q c^y* 1\T a yy\ & 




DETECT_MINAREA 


10 


DETECT_THRESH 


1.1 


DEBLEND_NTHRESH 


64 


DEBLEND_MINCONT 


0.005 


PH0T_AUT0P ARAMS 


2.5, 3.5 


PHOT_FLUXFRAC 


0.5 


BACK_SIZE 


64 


BACK_FILTERSIZE 


3 


BACKPH0T0_TYPE 


GLOBAL 



native scale of 0.05"pixel _1 to one of 0.03"pixel _1 , using 
the "Square" interpolation kernel. This is done primar- 
ily to better sample the ACS instrumental point-spread 
function (PSF). Lastly, the resampled image is registered 
to a common World Coordi nate System using positions 
from the USNO-B1 catalog (jMonet et al.ll2003l) . 

After processing the images, we create initial photo- 
metric galaxy catalogs f or each field using SExtractor 
(jBertin k, Arnouti 119961 ). As a goal of this study is to 
compare our results as closely as possible to those in G07, 
we use their SExtractor parameter set for our analysis. 
The parameters can be seen in Table [2] These parame- 
ters are selected to optimize detection of a suitably high 
number of small, dim objects, while also deblending close 
neighbors. This leads to a large number of spurious de- 
tections but the "false positives" are rejected during the 
analysis by applying a series of cuts to the data. This 
procedure is discussed more fully in the appendix. 

2.2. Spectroscopic Data 

While the determination of the weak lensing signal 
uses a generalized description of the redshift distribution 
of the background objects, the additional inclusion of 
the strong lensing information in the analysis (§3.3) is 
only possible if the redshifts of both the lens and the 
multiply-imaged background have been measured. Thus, 
to improve the utility of the lens sample, we obtained 
spectroscopic data on some of the best lens candidates 
from the COSMOS sample u sing the Low-Re solution 
Imaging Spectrograph (LRIS; lOke et all 119951) on the 
Keck I Telescope. The observations were conducted on 
UT 2009 Feb 22 and 23 in good conditions, with seeing 
ranging from 0'.'8 to 1". Spectra were obtained for 
the J095930.93+023427.7, 0038+4133, 0050+4901, 
0056+1226, 0211+1139, J100140. 12+020040.9, 
0216+2955, and 0254+1430 lens systems, where 
5930+3427, 0056+1226, J100140. 12+020040.9, 
0211+1139, and 0216+2955 were observed through 
slitmasks and the remaining systems were observed 
using a 1" longslit. The exposure times were set by the 
F814W magnitude of the lensing galaxy and ranged from 
1200 to 7200 s. The lensed source galaxies were much 
fainter than the lenses and, thus, the chosen exposure 
times were not sufficient to measure redshifts unless 
strong emission features or breaks in the spectrum were 
observed. In cases where we could not determine a 
redshift for the source galaxy, strong lensing constraints 
could not be obtained. However, these lenses are still 
useful for, and are included in, the weak lensing analysis. 
The data were reduced using custom Python scripts 



that performed the flat-field corrections, wavelength 
calibrations, rectifications, and extractions of the spec- 
tra. The extracted spectra were examined for multiple 
emission and/or absorption features. We were able to 
determine lens redshifts for six of the seven systems 
that were targeted (the slit for 0056+1226 was placed 
only over the lensed source because the lens redshift was 
previously known), but only two source redshifts. Thus, 
the new redshifts for COSMOS systems resulting from 
this work are: J095930.93+023427.7 (zi cns = 0.892), 
0038+4133 (z lcns = 0.738), 0056+1226 (z src = 0.808, 
based mostly on the 4000A break), J100140. 12+020040.9 
(z lons = 0.879), 0211+1139 (zi cns = 0.92, based mostly 
on the 4000A break), 0216+2955 (z lcns = 0.608), and 
0254+1430 (zi cns = 0.417, z SIC = 0.779). 

3. DISENTANGLING MASS AND LIGHT 

After reducing the data, we perform a full weak lens- 
ing analysis on the galaxy fields to infer the average mass 
properties of our lenses. The interested reader can find 
the full details of this analysis in the appendix. To sum- 
marize however, we measure the shapes of all background 
galaxies, apply a series of data cuts to remove contami- 
nants, and employ a system of tests to minimize system- 
atic uncertainty. We then convert these galaxy shape 
measurements into a measure of mass overdensity as a 
function of radius (AE(i?)), which we can compare with 
specific models to infer other properties of the lenses. 
The AE(i?) profile reflects the total mass distribution 
of the lens galaxy. We can investigate the dark mat- 
ter contribution by jointly modeling the luminous and 
dark components, constrained by the lens galaxy surface 
brightness profiles. To disentangle the luminous stellar 
mass profile from its surrounding dark matter halo, we 
employ a simple two-component mass model, whereby 
the stellar mass profile traces the light profile, and the 
dark profile exists as a single dark matter halo, centered 
on the galaxy itself (a one-halo central term in the frame- 
work of the "halo model" ) . This simplifying assumption 
explicitly excludes any mass contribution from an under- 
lying group/cluster halo, or nearby satellite galaxies (the 
two- halo term) , but we feel it is a well warranted assump- 
tion, because only a small fraction of massive elliptical 
galaxies (< 15% ) are located away from the center of 
their host halos (iMandclb aum et"aLll2006l) . 

3.1. Luminous Component 

We be gin our process of mo deling stellar light by using 
GA LFIlPl (iPeng et al.ll2002l) to fit a de Vaucouleurs pro- 
file (|de Vaucouleurslll948t) to the F814W images of each 
strong lensing galaxy, carefully masking out any contam- 
inating light from nearby satellites or (more importantly) 
the strongly lensed images. In most cases, this masking 
is achieved by running GALFIT on a small cutout of the 
galaxy (typically a square 2.5 times the galaxy's FWHM 
parameter - as determined by SExtractor - on a side). 
Nearly all of the light from our moderate-redshift lenses is 
contained well within the lensing Einstein radius, mak- 
ing it possible to exclude all of the lensed background 
structure without excluding significant portions of the 
light profile of the foreground galaxy. In the few cases 

12 http: / / users.ociw.edu / peng/work/galfit / galfit.html 
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Figure 1. Example of our stellar light profile model using GAL- 
FIT. Left: original data (COSMOS 0038+4133), showing the lens- 
ing galaxy and lensed Einstein ring. Middle: model of the lens 
galaxy. Right: residual image. Aside from a slight undcrsubtrac- 
tion in the core, the lens galaxy is well fit by the model. 

where this cutout method is not feasible (e.g. compact 4- 
image quasar lenses such as SDSS1138+0314 or HE1113- 
0641) we generate custom GALFIT masks to block out 
the remaining contaminating light. An example of the 
GALFIT modeling can be seen in Figured] 

By explicitly fitting with a de Vaucouleurs profile, we 
obtain two parameters for each galaxy that can be used in 
the mass profile modeling: an effective radius (R e ) and 
an F814W (Fband) magnitude (mi). We first convert 
the I-band apparent magnitudes into F555W (V-band) 
absolute magnitudes (My) by calculating K-corrections 
using the elliptical ga laxy template from the Kinney- 
Calzetti spectral atlas (|Kinnev et al.lll996|). and correct - 
ing for galactic extinction using the lSchleeel eTaTl (fi"998ft 
dust maps. This is done so that we may more readily 
compare our results to both G07 and previous mass-to- 
light (M/L) studies. To homogenize the sample, we pas- 
sively evolve all V-band lu minosities to a fiduc ial redshift 
z = 0.6, using the relation (jGrillo et al.ll2009D : 



d\ogL\ 

dz 



0.6 ±0.05 



(1) 



We also convert the total luminosity into effective sur- 
face brightness (I e ) using the equation: 



I, 



L v K 2 



2irR 2 n e K T(2n) 



(2) 



where n and k arc the general Scrsic parameters (for the 
de Vaucouleurs profile, n = 4, k = 7.67). The relevant 
parameters obtained through the photometric modeling 
are in Table [3] 

We model the stellar contribution to the mass profile 
as: 



S* = Tyl e exp 




- 1 



(3) 



where £* is the stellar mass surface density, and Ty 
is the rest-frame V-band stellar mass to light ratio 
(M*/Ly). Since both I e and R e are fixed by the GAL- 
FIT model for any given galaxy, this stellar mass profile 
has the benefit of being described by only a single free 
parameter, Ty. 

3.2. Dark Component 

For the dark matter, we assume a f unctional form 
of th e Navarro, Frenk and White (NFW; INavarro et al.l 
119971) profile. The projected, two dimensional NFW mass 
density profile (Enfw) can be described generally by: 

Tj NFW cx R s S c p c (4) 

(e.g., iWright fc Brainerdll2000l ). where R s is the NFW 
"scale radius", p c is the critical density of the universe, 



and S c = (A/3)c 3 /[ln(l + c) - j^] is a function of the 
concentration parameter (c = R&/R s ). i?A is the ra- 
dius at which the total mass density is A times p c . 
For each model, the value of A is determined solely by 
the redshift of the lens galaxy using the prescription of 
IBrvan fc Norman] (|1998l ) , which for our assumed cosmol- 
ogy allows i?A to be considered the virial radius of the 
system. For galaxies at z ~ 0.6, A is ~ 140, and all 
galaxies in our lens sample have overdensity values be- 
tween A ~ 125 and A ~ 160. 

By defining the virial radius, we are also able to fur- 
ther constrain the concentration parameter by assuming 
a functional form given by: 



1 



.12 x IQ^h-^M, 



© 



-0.14 



(5) 



as found in num e rical dark mat t er simulations 
(iBullock et al.1 I200H lEke et al.l I200TI : iHoekstra etafl 
I2005D . and where M vir is the total virial mass. While 
we assume no intrinsic scatter in the mass-concentration 
relation for our initial analysis, we do consider the effects 
of scatter, as well as other systematic effects that can im- 
plicitly alter this relationship, in §5.3. Rather than at- 
tempt to constrain the virial mass for each galaxy in our 
sample individually, we instead choose to parameterize 
these masses as a function of V-band luminosity (Ly). 
Specifically, we choose the form M V i r = TyLy, where Ty 
is defined to be the virial mass to light ratio. This then 
allows us to define the dark matter mass profile in terms 
of: 

SNFw(i?) = f(R, A, c) = f(R, z, M mr ) = f(R, z, L v , t v ) 

(6) 

As in the case of the stellar mass profile, we take Ly from 
the GALFIT models and, since the redshift is already 
known for each lens, we once again are left with a model 
with a single free parameter, this time Ty. 

3.3. Model Fitting 
3.3.1. Mass Overdensity Models 

We fit the models to the observed AE profile, de- 
termining the best-fit M/L ratios associated with the 
lens galaxies through a x 2 minimization procedure. We 
compare the observed weak lensing profile to the two- 
component models and average over contributions from 
individual models in order to mimic the procedure used 
in the weak lensing analysis. The merit function for our 
X 2 procedure is: 



Xwl 



'lte){ as,R - 



N,, 



J2 [AE*, j ( J R i ,Ty) + AS NFW> j(i?i,Ty)] i (7) 



3=1 



where AS(i?j) is the observed mass overdensity at ra- 
dius R4 (Table 0]), ctas.i is the error associated with that 
measurement, and AS„j and AS^fw j are the stellar 
and dark matter mass overdensity models evaluated at 
Ri (using the parameters of lens j), constructed accord- 
ing to Equation (|A12[) . This optimization to the weak 
lensing data set is performed for the full sample of lenses. 
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Table 3 

Photometric Lens Data. 



Lens Name 






mj 


K-corr 


E(B-V) 


My - 2.51ogi (fe) 


Ie 


(1 


PA 




(") 


I") 










(Lp. kpc -2 ) 






B0218+357 


0.17 


0.37 


19.89 


-0.19 


0.0680 


-20.10 


123.62 


0.97 


66.02 


SDSS0246-0825 


0.6 


0.99 


19.81 


-0.10 


0.0266 


-21.78 


557.53 


0.53 


24.88 


CFRS03P1077 


1.05 


1.89 


20.30 


0.42 


0.0983 


-23.46 


54.24 


0.66 


32.40 


HE0435-1223 


1.21 


1.12 


19.17 


-0.49 


0.0590 


-21.65 


99.52 


0.79 


-83.16 


B0445+128 


0.68 


2.29 


20.09 


-0.38 


0.3837 


-21.35 


12.94 


0.71 


-54.21 


B0631+519 


0.58 


0.84 


20.30 


-0.29 


0.0890 


-21.65 


104.31 


0.81 


8.35 


J0816+5003 


2.5 


0.94 


18.54 


-0.51 


0.0466 


-22.94 


293.74 


0.86 


46.23 


B0850+054 


0.34 


0.25 


21.73 


-0.34 


0.0615 


-20.60 


478.79 


0.42 


55.28 


SDSS0903+5028 


1.5 


0.55 


19.34 


-0.58 


0.0245 


-20.93 


293.68 


0.93 


-69.39 


SDSS0924+0219 


0.88 


0.58 


19.38 


-0.58 


0.0551 


-20.90 


248.06 


0.91 


52.62 


J1004+1229 


0.77 


0.51 


21.80 


0.45 


0.0372 


-21.88 


165.29 


0.38 


64.87 


HE1113-0641 


0.44 


2.35 


18.16 


-0.05 


0.0386 


-21.72 


168.03 


0.76 


-59.19 


Q1131-1231 


1.9 


3.63 


16.78 


-0.68 


0.0352 


-22.71 


55.83 


0.78 


10.93 


SDSS1138+0314 


0.67 


0.42 


20.12 


-0.50 


0.0196 


-20.60 


286.33 


0.79 


7.79 


SDSS1155+6346 


0.98 


0.52 


17.60 


-0.80 


0.0141 


-20.48 


898.45 


0.59 


36.23 


SDSS1226-0006 


0.63 


0.46 


18.98 


-0.42 


0.0233 


-21.32 


230.87 


0.28 


-34.38 


B1608+656 


1.14 


0.89 


19.51 


-0.28 


0.0310 


-20.24 


26.09 


0.44 


-13.98 


WFI2033-4723 


1.17 


0.73 


20.25 


-0.22 


0.0461 


-21.78 


140.47 


0.82 


41.79 


COSMOS5857+5949 


2.15 


1.23 


19.63 


-0.58 


0.0193 


-20.67 


45.30 


0.60 


58.41 


COSMOS5914+1219 


1.86 


1.27 


22.11 


0.75 


0.0202 


-22.33 


29.20 


0.69 


-22.72 


COSMOS5921+0638 


0.8 


0.33 


20.62 


-0.39 


0.0205 


-20.77 


366.25 


0.84 


-62.64 


COSMOS5941+3628 


1.21 


0.83 


20.90 


0.32 


0.0191 


-22.43 


120.08 


0.90 


69.28 


COSMOS5947+4752 


2.55 


0.35 


19.96 


-0.64 


0.0209 


-19.97 


365.42 


0.89 


89.20 


COSMOS0012+2015 


0.9 


0.51 


19.42 


-0.59 


0.0185 


-20.84 


321.34 


0.64 


68.23 


COSMOS0013+2249 


1.65 


2.02 


18.33 


-0.63 


0.0178 


-21.69 


52.91 


0.82 


-54.76 


COSMOS0018+3845 


0.4 


0.32 


23.13 


-0.13 


0.0187 


-19.20 


61.48 


0.68 


55.42 


COSMOS0038+4133 


0.74 


1.11 


20.36 


0.29 


0.0186 


-22.31 


82.28 


0.74 


88.38 


COSMOS0047+5023 


0.7 


1.33 


20.16 


0.19 


0.0182 


-23.06 


84.67 


0.80 


-56.86 


COSMOS0049+5128 


2.22 


0.31 


20.09 


-0.64 


0.0182 


-19.81 


428.99 


0.75 


22.46 


COSMOS0050+4901 


1.9 


0.74 


21.21 


0.48 


0.0189 


-22.53 


140.32 


0.73 


-65.24 


COSMOS0056+1226 


1.2 


0.76 


18.92 


-0.62 


0.0164 


-21.16 


210.72 


0.90 


-61.32 


COSMOS0124+5121 


0.84 


0.27 


22.33 


0.17 


0.0188 


-20.79 


271.75 


0.65 


38.48 


COSMOS0211+1139 


3.2 


1.66 


20.50 


0.32 


0.0168 


-23.05 


48.23 


0.61 


-83.30 


COSMOS0216+2955 


1.96 


0.99 


19.98 


-0.20 


0.0185 


-21.70 


84.55 


0.81 


62.14 


COSMOS0227+0451 


1.62 


0.71 


21.37 


0.29 


0.0176 


-22.05 


112.95 


0.60 


-10.32 


COSMOS0245+1430 


1.54 


1.70 


18.63 


-0.51 


0.0191 


-22.26 


72.62 


0.59 


31.75 


J095930.93+023427.7 


0.89 


0.95 


21.26 


0.28 


0.0191 


-22.07 


63.99 


0.69 


-2.33 


J100140.12+020040.9 


0.79 


0.29 


21.71 


0.26 


0.0180 


-21.58 


446.53 


0.83 


-71.80 


"Anchor" 


1.1 


0.49 


19.88 


-0.49 


0.0103 


-20.93 


285.63 


0.92 


-56.16 


"Cross" 


1.22 


0.90 


20.29 


0.09 


0.0085 


-22.60 


137.47 


0.79 


78.92 


"Dewdrop" 


0.76 


0.50 


19.96 


-0.35 


0.0089 


-21.59 


314.46 


0.85 


79.14 



Galaxy parameters obtained through GALFIT photometry fitting. Re%t\ is the Einstein radius of the lens (determined outside of GALFIT), and 
Re is the effective (half-light) radius of the de Vaucouleurs profile, mj is the modeled F814W apparent magnitude of the galaxy, K-corr represents 
the F814W to F555W K-corrcction, E(B-V) is the galactic extinction correction term, My is the absolute F555W magnitude (scaled by the Hubble 
parameter), and I e is the effective surface brightness of the galaxy's dc Vaucouleurs profile, as determined from My. Finally, q represents the axis 
ratio of the galaxy, and PA is its position angle, defined to be north through cast. While these final two parameters are reported as part of the 
GALFIT modeling, they are not used in the mass modeling of §3. 



To take full advantage of our data set, we also incor- 
porate the strong lensing information. We note that, for 
all strong lenses, the average value of the mass surface 
density within the Einstein radius is equivalent to the 
critical lensing density, £ cr it- This can be proven by not- 
ing that for a circularly symmetric lens, the tangential 
reduced shear profile (gt) can be written in the form: 



9t{R) 



Y,(<R)-Y>(R) 

S C rit — 



(8) 



is given by: 



Xsl 



« Ions, 2 

E 



E*(i?Ein,i) — £NFw(-REin,i) 



(9) 

where S C rit,i is the actual critical density of lens i, 
cScrit,! is the uncertainty on that value, £*(i?Ein,i) and 
^NFw(-REin,i) are the stellar and dark matter average 
mass density models of lens z, evaluated at the Einstein 
radius of that specific lens, and Ni ons ZB is the number of 
lenses for which the source rcdshift has been measured. 



Because gt is, by definition, equal to 1 at i?Ein it is a 
simple matter to show that £(< i?Ein) = ^crit- 

Thus, for all lenses for which we know both the lens 
and source galaxy redshift, we can accurately determine 
the value of the critical density of the system and further 
constrain the models. The strong lensing merit function 



3.3.2. Reduced Shear Models 

In addition to the mass overdensity model fit described 
in §3.3.1, we also fit models directly to our reduced shear 
data. Although shear is correlated with mass overden- 
sity, we do expect there to be a slight discrepancy be- 
tween the two methods: by creating models to fit reduced 
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shear directly, we no longer need to rely on the basic as- 
sumption of the weak lensing regime that convergence 
is always small (i.e., k <C 1). Indeed, this assumption 
will break down as we approach the Einstein radius, and 
to fully incorporate these smaller radii into the model, 
we should account for the non-linear response of lensing 
in this regime. To that end, we create two new merit 
functions that parallel the Xwi an d Xsi used in the mass 
overdensity fit. 
For weak lensing data, we use the function: 



Xwi, shear 



E 



i 



9t(Ri) 



Ni, 



St,', 

Ni cns ~. 2 



, N lcns 



3=1 



where gt(Ri) is the observed reduced shear signal in a 
given radial bin (Tabled, Cg t ,i is the measured uncer- 
tainty, and 3t mod .j is the model combined (de Vaucouleurs 
+ NFW) reduced shear profile for lens j, given by: 



9t n 



J NFW 



) — (S* + Enfw) 



E C rit 



J NFWi 



(11) 



since individual components of reduced shear do not add 
linearly. Once again, xt\ shear ^ s calculated for the full 
sample. 

For lenses where both Z\ and z s are known, we create 
a strong lensing merit function, of the form: 



Xsi. shear 



= E 



REin,i — REm{gt mod ,i) 



(12) 



where -REin,i is the observed Einstein radius of lens i, 
o~R ELn ,i i s ^s measured uncertainty, and i?Ein(st mod ,i) is 
the model Einstein radius - determined by inverting 
Equation (|lll) to find the radius where the model reduced 
shear profile is equal to 1. 

4. RESULTS 

4.1. Total Mass Profile 

The observed mass overdensity profile can be seen in 
Figure [2] We can see a significant signal in the AS (E- 
modc) profile, while the AS X (B-mode) profile is consis- 
tent with zero. Measured values of both AS and AS X 
are reported in Table |H along with their associated shear 
values. To measure the slope of the density profile, we 
fit both a generic power law model (AS oc R 7 ) and an 
isothermal model (AS oc Rr 1 ) to the data (Figure^]). We 
find that the best-fit power-law model (7 = — 0.86±0.16) 
does not deviate from the best-fit isothermal model, pa- 
rameterized by a velocity dispersion of o~ v — (260 ± 20) 
km s" 1 , by more than 1-er. Given our uncertainties, both 
models are consistent with the data, allowing us to con- 
clude that the average mass overdensity profile of our 
sample is approximately isothermal. 

4.2. Luminous and Dark Matter Profiles 

We minimize the total merit function for each of the 
two approaches: Xtot = Xwi + Xsi f° r the mass overden- 
sity fit, and Xtot.shear = Xwhshear + Xs^shear for the re- 
duced shear fit. The best- fit M/L parameters obtained 



with both approaches can be seen in Figure [3J The re- 
sults from the model fitting procedures disagree by less 
than one standard deviation. Thus, we feel we can use 
the results from the mass overdensity model fit without 
worrying about inaccuracies in the profile due to the sim- 
plifying assumptions made about shear. 

The mass overdensity profiles generated from our best- 
fit parameters are shown in Figure [4J plotted along with 
the total mass profile from our weak lensing analysis. 
The figure shows the contributions of the stellar (blue 
dashed line) and dark matter (red dash-dotted line) mass 
profiles to the best-fit total mass profile (solid black line). 
The best fit to our data has xl = 0.98, indicating a good 
agreement between our model and the observed overden- 
sity profile. As a comparison, the pure power-law fit in 
§4.1 has xt — 0.81, meaning that, given our errors, nei- 
ther the the 1- nor the 2-component model fit is strongly 
preferred over the other. 

We find a virial M/L ratio t v = 300 ±90 h M e /L e 
and a stellar M/L ratio Ty = 4.0 ±0.6 h Mq/Lq, for an 
overall virial to stellar mass ratio of M V i r /M* = 75 ± 25. 
Given the mean V-band luminosity of our sample, Ly = 
4.5 x 10 10 h~ 2 Lq, this corresponds to average stellar 
and virial masses of M* = (1.8 ±0.3) x 10 11 h^ 1 M Q and 
M vir = (1.4 ±0.4) x 10 13 h" 1 M Q respectively. Further- 
more, the average virial mass we find corresponds to a 
typical virial radius of 332 hT x kpc and a typical NFW 
scale radius of 62 h~ x kpc. 

Combining our mass ratio data with the assumed cos- 
mology, we are also able to estimate the stellar baryon 
fraction (/*) present in our lens sample, according to the 
relation: 

/, - (13) 

where the quantity represents the total baryon to 
dark matter ratio in the galaxy, and is given by the 
globa l value of ^ = 0.176 ± 0.013 (jSpergel et al.l 
120071) . For our sample, we find a stellar baryon frac- 
tion of /* = 0.075 ± 0.030, wh i ch is in excellent agree- 
ment with IMandelbaum et aLl ([2006D . who find /* = 
0.055!qq1o l° r their "sm6" bin, consisting of early- 
type galaxies having M* = 2.13 x 10 11 h~ x M© and 
1.58 x 10 13 Mq, and in good agreement with 
who find /* = 0.10 ± 0.03. The 
galaxy sample is somewhat differ- 
ent from our own, in that it consists of both early- and 
late- type galaxies. The sample is dominated (~ 75% ) 
by early types, though, and the authors note that adopt- 
ing a more restrictive selection criterion of Sersic index 
n > 2.5 does not change t heir results. Th e aver age stel- 
lar and virial masses of the lHevmans et aLl (|2006l ) sample 
(~ 7x 10 10 h- 1 M Q and - 3 x 10 12 h' 1 M Q , respec- 
tively) are an order of magnitude smaller than our own 
sample, making it difficult to accuratel y compare results. 
However, as the iHevmans et aLl (l2006fl sam ple has simi- 
lar masses to the IMandelbaum et aLl (|2006h "sm4" early- 
type sample, scaling this value by /*, S m6//*,sm4 = 0.33 
gives a result of /* = 0.03 ± 0.02, which is still within 
1.25 — a of our value. 

Finally, our val ue is in marginal agreement with 
iFukugita et aLl (|1998f ). who estimate a universal stellar 
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□ B-mode Overdensity Profile 
• Mass Overdensity Profile 

- SIS Profile a = 260 km/s 

- Power Law Profile 7 = —0.86 




10 

Radius [l&c/h] 



1 Mass Overdensity Profile 
SIS Profile a = 260 km/s 
Power Law Profile 7 = -O.t 
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Radius [kpc/h] 



Figure 2. Mass overdensity profile. The best-fit generic power law model (red dashed line) is consistent with the data, and does not 
significantly deviate from the best-fit isothermal profile (green line). 



Table 4 

Measured Excess Surface Density and Tangential Shear. 



Galaxy Counts 


Outer Bin Boundaries 21 
(h- 1 kpc) 


AS 




0"AS 


at 


9x 


°"9 


4 


11 


1213 


13 


580 


0.231 


-0.010 


0.101 


23 


19 


755 


110 


228 


0.124 


0.023 


0.049 


73 


33 


108 


-73 


130 


0.024 


-0.010 


0.027 


233 


57 


140 


-102 


69 


0.029 


-0.020 


0.015 


656 


100 


90 


66 


43 


0.017 


0.014 


0.009 


1851 


176 


56 


2 


25 


0.012 


0.000 


0.005 


4380 


308 


38 


14 


16 


0.007 


0.004 


0.004 


8287 


541 


33 


-6 


12 


0.008 


-0.001 


0.003 


6997 


949 


23 


8 


15 


0.005 


0.001 


0.003 



A quantitative description of the mass overdensity profile (Figurc[2j for our lens sample. Here, Galaxy Counts refers to the number of galaxies 
in each radial bin (given by projected radius). The outer radius of each bin is shown in the radial bin boundaries column. The next three columns 
show respectively the excess surface mass overdensity, "B-Mode" overdensity, and uncertainty on each mass measurement in a given bin. The final 
three columns show the associated shear quantities. 
a The inner radius for the innermost bin is 6 h~ 1 kpc. 
800 r 



700 




Figure 3. Confidence contours for the best-fit M/L ratios for 
mass overdensity models (red circle) and the associated 68%, 95%, 
and 99.7% confidence regions. The best-fit value for the reduced 
shear analysis is shown by the blue square. The results of the fits 
do not significantly disagree with one another. 

baryon fraction of /* ~ 0.12 for spheroids, using data 
acquired from the local universe. 

5. DISCUSSION 
5.1. Mass Profiles 



t v = 300 hMJL® 
~£ v = 4.00 hM r JL e 



' Observed Data 
Model deVauc Profile 
Model NFW Profile 
Total Mass Profile 




Figure 4. Stellar and dark matter profiles generated by our best- 
fit mass-to-light ratios. The blue dashed line represents the stellar 
de Vaucouleurs density profile, while the red dash-dotted line rep- 
resents the dark matter NFW profile. 

Figure [2] shows that the best-fit power law profile to 
the total mass overdensity data between physical scales 
of -10 and -1000 h" 1 kpc has 7 = -0.86 ± 0.16, where 
AS(i2) oc R 1 . This is consistent with an isothermal 
(7=1) profile and, thus, is similar to the results of 
G07, who found that the low-redshift SLACS lens sample 
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Figure 5. Top: Evolution of the ratio of dark to luminous matter 
between SLACS and our lens sample, after splitting our sample into 
two separate redshift bins. The best-fit slope of the data is shown 
as a solid line, with a = 51±36 and = 0.9±1.8. Bottom: Redshift 
histograms for each data point, showing the redshift distribution 
of lens galaxies within each bin. 



Table 5 

Mass Ratio Evolution. 



Sample 


z 




M vir /M* 


a M vir /M» 






(10 10 h~ 2 Lq) 




SLACS 


0.23 


5.8 


54 


28 


zl 


0.42 


1.2 


85 


33 


z2 


0.83 


■1.8 


80 


46 



Data values associated with Figurc[5] The z and L v values for each 
data point arc, respectively, the average redshift and V-band luminosity 
of all lenses included in the bin. 

also displays a characteristically isothermal mass profile 
between ~1 and ~300 hr l kpc. However, a closer in- 
spection of the right panel of Figure [D reveals that we 
may be seeing some flattening in the profile at large 
radii. These points are located at distances outside of 
the average virial radius of our lens sample. Background 
galaxies located at these radii could thus be affected by 
the dark matter halos of other massive galaxies nearby 
the lens, which would explain the larger-than-expected 
shear signal in the profile. Alternatively, since many of 
the galaxies in these bins lie near the edge of the ACS 
field, we could be seeing excess "shear" due to CCD edge 
effects. Rerunning the profile analysis after excluding 
the outermost two data bins, we find a best-fit profile 
7 = —1.06 ± 0.09 that is much more consistent with an 
isothermal distribution. 

Even with the possibility of flattening at large radii, 
though, combining the results of G07 with our results 
suggests that the average mass overdensity profile of 
strong gravitational lenses is consistent with an isother- 
mal model over three orders of magnitude in physical ra- 
dius. This result is in good agreement with previous mass 
profile studies of non-lensing ellipticals (|Sheldon et al.l 
[200l iMandelbaum et al.ll2006l ) between - 50 - 300 h' 1 
kpc, which would suggest that there is little difference 
between the mass profiles of field strong lenses and their 
non-lensing counterparts at moderate radii. 

Comparing the mean velocity dispersion of our lens 
sample (a% = 260 ± 20 km s" 1 ) to that of the G07 SLACS 
sample (a^ ~ 248 km s _1 ), and noting the broad similar- 



ities between the slopes of their mass profiles, it would 
appear that there is little difference between the overall 
mass properties of strong lenses at z = 0.6 and the prop- 
erties of those at z = 0.2, although large uncertainties 
on both profiles could easily hide any real evolutionary 
results. Additionally, the best- fit M/L ratios of our sam- 
ple (r y = 300±90 h M & /L Q , T v = 4.0±0.6 h M e /L Q ) 
are not significantly different from those of G07 (ry = 
246 ± 100 h Mq/L q , Ty = 4.48 ± 0.46 h M Q /L Q ). 

Considering the hierarchical merging model of galaxy 
formation, which posits that the majority of the mass 
assembly of massive red ellipticals (the dominant mor- 
phology of lensing galaxies) takes place at redshifts z > 
2, it is not surprising that the differences in total mass 
between these two samples is not dramatic. More sur- 
prising however is the fact that the individual M/L ra- 
tios do not change, as passive luminosity evolution could 
significantly affect the total brightness of a galaxy over 
this redshift range. We investigate the effects of purely 
passive evolution on the stellar M/L component in two 
ways, back-evolving the G07 Ty res ult from z = 0.2 to 
z = 0.6 using b oth the methods of iTreu et all (|2001l ) 
(method 1) and Ivan Dokkum fc Franxl ^OOiF Tmethod 
2). After applying these relations, we find a theoretical 
M/L ratio of Tv. z =o.6 = 3.1 ± 0.3 using method 1 and 
T v,z=o.6 = 3.4 ± 0.3 using method 2. Both results agree 
with our measured Ty within the errors. From this, 
we conclude that the passive evolution and no-evolution 
models are both plausible descriptions of stellar M/L ra- 
tio evolution between ~z ~ 0.2 and z ~ 0.6, given the cur- 
rent size of our errors. Future studies with larger sample 
sizes should be able to place tighter constraints on these 
parameters, making it easier to distinguish between the 
possible outcomes. 

5.2. Mass Ratio Evolution 

We also search for evidence of the evolution of mass 
ratio, looking for any significant changes in the ratio of 
M V i r /M*, which could suggest evidence of a merger or 
satellite accretion event (jConrov et al.ll2007h . We do this 
by comparing our results to those of G07, after splitting 
our sample into two redshift bins: subsample zl, where 
z < 0.6 and subsample z2, where z > 0.6. We fit a slope 
to the data of the functional form: 



-(z) =a(l + zf 



(14) 



The results of the fit can be seen in Figure[Sl and the data 
are presented in Table [5] We find a best-fit a = 51 ± 36 
and a best-fit = 0.9 ±1.8, which is consistent with a no- 
evolution model (0 = 0). This res ult is in agreement with 
the work of lHevmans et all (|2006D . who use weak- lensing 
alone to probe the evolution of a sample of predominantly 
(~ 75%) early- type galaxies between z = 0.2 and z = 0.8. 
After explicitly assuming a linear growth factor, they find 
that M vir /M*(z) = (3 4±ll) + (31±35 ) z. In using weak- 
lensing alone however. iHevmans et al.l (|2006|) are unable 
to constrain stellar mass directly, since weak-lensing is 
not able to probe mass at small radii where stellar mat- 
ter dominates. Instead, they make an assumption about 
the nature of stellar mass, using a Kroupa stellar IMF 
(|Kroupa et al.l [l993h to scale the average halo mass for 
their galaxy sample. This is similar to techniques used 
by other weak-lensing only mass ratio studies, such as 
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Figure 6. Relationship between inner dark matter slope 7 and 
stellar M/L ratio Ty, showing a degeneracy between the param- 
eters and demonstrating the fact that increasing the dark matter 
fraction in the core necessarily decreases the stellar mass fraction, 
assuming a constant total mass. The data point represents our 
best fit Ty value assuming an NFW (7 = 1) profile and the errors 
represent the 1-u confidence interval of that fit. The dashed lines 
show the limits on Ty when 7 varies by 20%. 



Mandelbaum et all (120061) who used a Kroupa IMF and 
Hoekstra et al.l (|2005[ ) who used a scaled Salpeter IMF 
(lSalpetenll955fi . 

Unlike these previous studies, our stellar mass values 
are obtained purely from a luminosity-scaled fit to the 
mass data, without the need to invoke a specific stellar 
IMF. Indeed, this is the first such study to attempt to 
quantify the evolution of the stellar mass ratio without 
making any assumptions about the nature of the stellar 
IMF, and the fact that we are able to obtain a consis- 
tent result, even with such a small sample of lenses is 
promising. It is likely that future studies with larger 
lens samples will be able to place even better constraints 
on stellar mass evolution, without having to worry about 
possible systematic errors associated with the IMF, fo- 
cusing instead on only the lensing-inferred mass and total 
luminosity, which are both directly determined from the 
data. 

5.3. Mass to Light Ratio Systematics 

In determining the best-fit stellar and virial M/L ratios 
for our lens sample, we determine statistical uncertainties 
from the model-fitting procedure that are on the order 
of ~ 20 — 30%. In this section, we investigate the impact 
of systematic uncertainties. Specifically, we look at two 
such effects: varying the inner slope of the dark matter 
profile, and modifying the dark matter concentration- 
mass relation (CMR). 

5.3.1. Dark Matter Inner Slope 

When discussing the disentanglement of stellar and 
dark matter, we have, to this point, specifically assumed 
a standard NFW profile for the dark matter. This 
has been done not only to take advantage of the ana- 
lytic form of the dark ma tter shear profile described in 
iWright fc Brainerdl (|2000t ). but also to compare our re- 
sults to the numerous previous studies th at have assumed 
NFW dark matter profile s as w el l (e.g. iHoekstra et alj 
120051 iMandelbaum et al.1 [20061: iHevmans et al.1 120061: 
G07). Recently however, studies have shown that lens- 
ing data actually favor a dark matter profile with a 



slightly modified inner slope (7), often motivated by 
some physical process such as adiabatic cont raction (AC; 
IBlumenthal et al.|[l98l iGnedin et all 120041) that affects 
the re gions of a galaxy where R < R*j Jiang fc Kochanekl 
120071 lGaoiraLll200a fSchulz et al.l l2009f >. Because this 
modification only alters the distribution of the dark mat- 
ter on small (< 50 h~ x kpc) scales for galaxy-scale 
masses, the total dark matter mass and virial M/L ra- 
tio derived from these new profiles remain statistically 
consistent with the original NFW model. However, this 
is not the case for the derived stellar M/L ratio (T), as 
even a change in 7 on the order of a few percent can sig- 
nificantly impact T, and hence the inferred stellar mass. 

To gauge the systematic impact of altering 7 on our 
data, we rerun the two-component stellar+dark matter 
analysis described in §3, replacing the standard NFW 
profile with a generalized NFW (gNFW) profile (e.g. 
lKeetoril200ll) : 



p(r) = 



(r/r,)T(l + r/r a ) 



3-7 



(15) 



that is similarly projected into 2-dimensions fsee lKeetonl 
120011 for details). While this prescription is somewhat 
different than applying an AC model to a galaxy (AC 
also modifies the concentration parameter, whereas al- 
tering the gNFW slope does no t), the effect on the halo's 
inner slope can be comparable: IGnedin et al.l (|2004l ) have 
demonstrated that, assuming an initial NFW dark mat- 
ter halo with a distribution of baryons that condenses to 
form an elliptical galaxy, the final inner slope of the dark 
matter halo ranges between 7 ~ 1.1 and 7 ~ 1.5. 

To be conservative, we vary 7 between 0.1 and 1.5, 
measuring the best- fit V-band stellar M/L ratio (Ty) 
for each case. The results can be seen in Figure [6] We 
find that the best-fit Ty is strongly dependent on 7, sug- 
gesting a degeneracy between these two parameters that 
affects stellar mass in the expected way: an increase in 
dark matter fraction in the core of the galaxy (character- 
ized by a steeper inner slope) results in a smaller fraction 
of stellar mass (characterized by a lower stellar M /L ra- 
tio). 

Assuming a fiducial uncertainty on 7 of 20%, the vari- 
ation in 7 corresponds to a range in stellar M/L ratio 
of 2.5 h (Mq/Lq) < Ty < 4.8 h (M®/L @ ). Adding 
this systematic uncertainty to our previous estimate of 
statistical error yi elds (Ty = 4.0 ± 6 h M Q /L Q ). 
As a comparison, Uiang fc Kochanekl (|2007f l measure a 
~30% decrease in stella r M/L ratio by including the 
IBlumenthal et al.l (|1986h AC model on a subsample of 
lower-redshift strong lenses, in good agreement with the 
results we see when we increase the inner slope of the 
NFW profile by -20%. 

5.3.2. Concentration-Mass Relation 

By using a CMR to reduce the number of free pa- 
rameters in our model fitting routine, we are explic- 
itly coupling the best-fit mass to the observed shape 
of the mass density profile. If instead wc were to 
make different assumptions about the CMR, this would 
lead to different best-fit mass values, which ultimately 
would yield systematically different stellar and virial 
M/L ratios. While we, like G07, have assumed the 
IBullock et ail (|200l CMR (Equation ©) in our anal- 
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ysis, recent work involving higher res olution numerical 
simulations (|Macci6 et al.ll2007l 12008ft have argued for 
an alternative CMR that is less sensitive to halo mass 
and evolves more slowly with redshift. In particular, us- 
ing the c osmological pa r amete rs from the WMAP 5-year 
data set, IMaccio et all (|2008l ) determine a CMR of the 
form: 



9.354 



M„ 



(H(z)/H ) 2 / 3 \l x lO 12 /i" 1 M 



-0.094 



(16) 



where H(z) is the Hubble parameter. Comparing these 
CMRs for galaxies with total mass of the order of our 
sampl e (~ 10 13 M Q ), we find that the IMaccio et al.l 
(|2008f) CMR yields a co ncentration paramet er that is 
~ 15% smaller than the iBullock et at] (|2001ft CMR at 
redshift z = 0, but (because of the slower evolution 
with redshift) actually becomes ~ 10% larger than the 
IBullock et ail (|2001l ) CMR at z = 0.6, the median red- 
shift of our sample. Re-running our initial de Vau- 
couleurs + NFW model fit with this new CMR, we find 
that the best fit ry and Tv parameters vary by approxi- 
mately 10% and 5% respectively, which is much less than 
the statistical uncertainties presented in §4. 

In addition to completely modifying the CMR, we 
also investigate systematics associated with the intrin- 
sic scatter of the rela tion itself. Including the 1-er errors, 
IBullock et all (l200l measure a ~ 40% variation in halo 
concentrations for their sample of ~ 10 13 M Q ha- 
los, which we incorporate into our analysis. Specifically, 
assuming an extreme case where we increase the concen- 
tration parameter of each galaxy in our sample by ~ 40% 
(and the alternative case where we decrease each halo by 
~ 40%), we measure variations in ry and Ty of 25% and 
15% respectively. Adding these uncertainties in quadra- 
ture with the uncertainties associated with our choice 
of CMR, we therefore include an additional ~ 15% and 
~ 25% uncertainty on our best-fit stellar and virial M/L 
ratios due to the systematics of the CMR. 

5.3.3. Propagation of Systematic Uncertainties 

Ultimately, by investigating modifications to both the 
inner slope and concentration parameter of the dark mat- 
ter halo, we believe that we are able to characterize the 
potential impact of important systematics on our results, 
enabling us to propagate these uncertainties to all other 
measured quantities that depend on them. In particu- 
lar, we find a total systematic uncertainty of ~ 30% on 
our best-fit Ty (dominated by the systematics of the 
dark matter inner slope) and a ~ 25% uncertainty on 
our best- fit Ty , which subsequently correspond to ~ 30% 
and ~ 25% uncertainties on our derived stellar and virial 
masses. Propagating these uncertainties to the other de- 
rived parameters, we measure a ~ 40% systematic un- 
certainty on both the mass ratio {M V i r / M*) , and stellar 
baryon fraction /». 

6. CONCLUSIONS AND FUTURE WORK 

Taking advantage of the large density of background 
galaxies obtained by relatively short (~ 1800 sec) ACS 
exposures, we present the results of a joint weak and 
strong lensing analysis of a sample of 41 massive ellipti- 
cal galaxies at moderate (z ~ 0.4 — 0.9) redshift, deter- 
mining the mass overdensity profile over nearly 3 orders 
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Figure 7. Comparison of total luminosity (extinction corrected 
and passively evolved to z = 0) vs. observed redshift between the 
deep-exposure SLACS lens sample and the lenses used in this work. 
The SLACS lenses are more luminous than our sample. 



of magnitude in physical radius. Using a 2-component 
de Vaucoulcurs + NFW profile model, we are able to 
determine stellar and virial M/L ratios for the sample, 
and disentangle the relative contributions of stellar and 
dark matter from the total mass budget. Furthermore, 
we compare all of our results to those obtained from a 
subset of the lower redshift SLACS sample (z ~ 0.2), 
placing constraints on the evolution of these mass prop- 
erties over cosmic time. Our results can be summarized 
as follows: 

• We present new redshift information for 
7 COSMOS lenses: we find lens red- 
shifts for J095930.93+023427.7, 0038+4133, 
J100140.12+020040.9, 0211+1139, and 
0216+2955, a source redshift for 0056+1226, 
and both lens and source redshifts for 0254+1430. 

• The mass overdensity profile of our sample has a 
best-fit power-law profile of 7 = —0.86 ± 0.16 and 
is consistent with an isothermal model between ~ 
10 h,- 1 kpc and - 1000 h,- 1 kpc. 

• The best-fit stellar and virial M/L ratios for 
our lens sample are given by Ty = 4.0 ± 
0.6 ti.l h Mq/L q and t v = 300+90+75 h Mq/Lq 
respectively. 

• Given our average sample luminosity of Ly — 
4.5 x 10 10 h~ 2 Lq, wc find an average virial mass 

(1.4 ± 0.4 ± 0.3) x 10 13 h' 1 Mq and 

(i.8±o.3 ±° : 5 8 )x 



of M vir 
an average stellar mass of 
10 11 h- 1 Mcr 



• We find an average mass ratio of M w ,>/M, = 75 ± 
25 ^27 1 which is consistent with the lower- redshift 
SLACS sample ratio of M vir /M* = 54 + 28. 

• Using our assumed cosmology, we convert mass 
fraction into a stellar baryon fraction, finding that 
/* = 0.075 ± 0.030 toMQ- This result indi- 
cates a low stellar formation efficiency in early- 
type galaxies, in agreement with previous resu lts 
(jMandelbaum et al.ll2006t iHevmans et al.ll2006l ). 
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• By comparing our mass fraction data to the lowcr- 
rcdshift SLACS sample, we are able to place con- 
straints on the evolution of the virial to stellar 
mass fraction of massive early- type galaxies. We 
find that the quantity evolves as M V i r (z) / M*(z) = 
(51±36)(l + z)(°- 9±1 - 8 ) over the last - 7 Gigayears. 
This is the first such study to place constraints on 
the evolution of M V i r /M t without invoking a spe- 
cific IMF. 

Overall, these results show that our method is a 
promising new technique: with only a small sample of 
lenses (41 systems), were are able to significantly mea- 
sure the average mass overdensity profile and mass-to- 
light ratios of the sample, allowing us to place constraints 
on the mass evolution of early-type galaxies from z ~ 0.8 
to z ~ 0.2. As always though, this work can be strength- 
ened by increasing the sample size. While galaxy-scale 
strong lenses are relatively rare today (~ 100 systems), 
future large-scale observational surveys such as LSST 
(~ 10000 new lenses ) and JDEM (~ 10000 new lenses; 
iMarshall et al.l 120051 ) are expected to dramatically in- 
crease the number of known strong lenses. Including 
these new lenses in future work will greatly reduce the 
size of our present statistical uncertainties. 

In the near future however, we plan to increase our 
lens sample by including other known deep space-based 
lens data from CASTLES that have been imaged us- 
ing the smaller HST Wide Field / Planetary Camera 2 
(WFPC2). Though the smaller area will limit our ability 
to probe out to the large radii made available with the 
ACS sample, the longer exposure times for the WFPC2 
data yield background densities that are comparable to 
those in the ACS sample, allowing us to place better 
constraints on the inner mass overdensity profile by in- 
corporating these galaxies into our current lens sample. 
We will also take advantage of the HAGGLeS stro ng-lens 
search of bright red galaxies ((Marshall et al. 2009) which 
is expected to explore the whole of the HST ACS archive, 
discovering ~ 10 lenses per square degree of coverage. 
We plan to further augment this study by combining our 
ACS sample with the full (~ 40) long-exposure SLACS 
sample, which will allow us to probe the differences in the 
mass properties of strong lenses across a wider variety of 
samples and categories, such as morphology, luminosity, 
and finer redshift slices. We are particularly interested 
in segregating the lens sample by total luminosity, since 
the SLACS lenses appear to be more luminous than their 
modcratc-rcdshift counterparts when passively evolved 
to z — (Figure [7J| , suggesting that the SLACS sample 
may have evolved from a slightly different population of 
galaxies. 

Of course, by increasing the sample size of lenses and 
thus reducing the statistical errors associated with these 
measurements, correcting systematic uncertainties will 
become much more important. In future studies we 
hope to improve our control of systematics in several 
ways. Foremost is the ability to measure accurate red- 
shifts for both the lens and source galaxies of strong lens 
systems, as these measurements are crucial for accurate 
lens modeling. While some future survey instruments 
should provide redshift information (cither spectroscop- 
ically or photometrically) as part of their normal oper- 
ation, it is still important to improve the redshift infor- 



mation for currently known lenses, and we plan to obtain 
more spectroscopic and photometric redshift data in the 
near future. Additional systematics controls could come 
in the form of improved profile modeling, relaxing the 
strict de Vaucoulers + NFW mass profiles in favor of a 
more general free-index Scrsic profile for the stellar mass 
and a freely- varying gNFW profile for the dark matter. 
Similarly, improved strong lens mass modeling should be 
utilized in the future, increasing our ability to distinguish 
between various mass profiles at very small scales, as well 
as improving the precision of our stellar mass fraction 
measurements. Finally, we could improve galaxy selec- 
tion and weak lensing measurements by using HST data 
obtained through other filters, allowing us to determine 
photometric rcdshifts for many of the background source 
galaxies. We note however that this should not be a 
dominant source of systematic error, as we have shown 
previously that altering the background galaxy redshift 
distribution does not significantly change the results of 
any of our fitted parameters. 
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APPENDIX 
WEAK LENSING ANALYSIS 
Object Detection and Shape Measurement 

This study utilizes the IMCATr 1 ^! software suite, which is based on the well-known algorithm of Kaiser, Squires, & 
Br oadhurst (1995, hereafter KSB), as well as a series of Perl and Python scripts which apply the PSF correction scheme 
of iSchrabback et al.l ([2007D and mass profile analysis of G07. Using both the catalogs and their associated science 
images, our pipeline first invokes the IMCAT object detection program hfindpeaks on the images and cross-correlates 
the results with the related SExtractor catalogs, in order to reduce the number of false detections. Objects detected 
in both SExtractor and hfindpeaks are convolved with a Gaussian smoothing kernel in order to regularize their shapes 
and sizes and also to reduce shape noise; the width of an object's smoothing kernel (r g ) is given by the object's half 
light radius as determined by SExtractor. Once an object has been smoothed, the local sky background is calculated 
and subtracted from its overall intensity. After ba ckg round subtraction, the object's magnitude is determined using a 
magnitude zero point taken from the ACS websitcO- 

Finally, the script determines a centroid for each object, and calculates the quadrupolc moments of intensity, given 
by 

Qij= f d 2 8W(8)9i6jI(8), i,j£x,y (Al) 

J r<4r g 

where the integral is taken from the centroid location to 4 x r g . Here, 1(9) represents the intensity of an object in a 
given pixel, and W(6) is a weight function, which for this study is a circular Gaussian with radius r g , matching the 
filter used for object detection. These moments of intensity are then used to calculate "ellipticity polarizations" : 

62 = <TW" (A3) 

^xx ~ ^yy 

which represent the components of the observed shape of the object. In addition, the "sm ear" and "shear" pola rizability 
pseudo-tensors ( y P sm and P sh respectively; see e.g. KSB, or the updated definitions in lHoekstra et al1l!998| ) are also 
calculated, which are used to correct shape measurement errors due to PSF distortion. 

Charge Transfer Inefficiency Correction 

Before we analyze the weak lensing signal, we must correct for any systematic distortions of the galaxy shapes. The 
first such systematic that we address is distortion due to the degradation of charge transfer efficiency (CTE) of the 
ACS CCDs. Defects that arise during the lifetime of the CCD can cause some charge to be delayed and transferred to 
other pixels, resulting in a trail of charge following an object. In the ACS, the readout direction is vertical. Thus, this 
trail will induce a spurious stretching in the — e\ direction. Since CTE degradation alters the shapes of galaxies in the 
same direction, the effect can mimic the gravitational shear signal caused by weak lensing. The spurious polarization 
induced by CTE degradation is typically on the order of 1%, and therefore it most strongly affects the outer regions of 
a galaxy-galaxy shear profile (where the gravitational shear drops below 1%). If left uncorrected, this spurious signal 
can alter the inferred total mass of the galaxy. 

The "strength" of an object's tail is affected by several factors, including object brightness, size, and location along 
the CCD readout path. This makes for an in herently non-l i near p roblem that cannot be corrected by a simple model. 
Instead, we use the empirical prescription of iRhodes et al.l (|2007|) devised for the COSMOS program. Following G07, 
we modify the original equation to be compatible to our data, although we use a slightly different prefactor which 

13 http://www.ifa.hawaii.edu/~kaiscr/imcat/ 14 // . . j /, . / / 

1 1 ! 1 11 http://www.stsci.edu/hst/acs/ 
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Figure 8. Average e\ component of ellipticity as a function of relative y pixel (j/), with y' = [y' = 1] representing the bottom [top] 
of the ACS field of view. The effect of CTE degradation is seen as an average nonzero ellipticity in the — ei direction, and is seen most 
prominently in our uncorrected ellipticity bins (red squares) near the region farthest away from the CCD readout registers (y' = 0.5), as 
we expect. Our correction for CTE degradation (blue circles) removes nearly all of the spurious stretching. 

accounts for combining the IMCAT de finition of SNR with the SExtractor measured half-light radius (G07 obtain 
both parameters from IMCAT, whereas [Rhodes et al.l (|2007t ) use SExtractor exclusively): 



eiCTE 



-6.5 x 10" 



2- y 



(SNR) 



-0.9 



x (MJD - 52333) ( 



n 



0.18' 



(A4) 



which modifies the e± component of each object in a catalog. In this equation, represents the half-light radius of 
the object (as determined by IMCAT), and y' is the normalized y-coordinate of the object (equal to at the bottom 
of an ACS field and 1 at the top), defined such that the correction for the CTE degradation is maximized near the 
ACS chip gap (y 1 ~ 0.5), where the pixels are the farthest away from the readout registers. The time-dependent factor 
takes into account the fact that the effects of CTE degradation have steadily worsened over the lifetime of the ACS. 

Figure [5] shows the effect of CTE correction on our data sample, and from the figure we can see that the average 
corrected galaxy shape (ei^orr = ei — ei.cTE) is no longer spatially dependent, and is very nearly zero throughout 
the entire field. In addition, our CTE correction scheme is able to correct specifically chosen subsets of the galaxy 
population as well: after separating the source galaxies into three distinct magnitude bins (m < 25; 25 < m < 26; 
m > 26), we find that the CTE-corrected galaxies in each bin show a mean shape ei jCorr = 0, again with no spatial 
dependence across the field. 

To test the systematics associated with our CTE correction scheme, we apply an additional 20% CTE shape correction 
to every background galaxy, regardless of size, brightness, or position on the detector. We perform a full weak 
lensing analysis and 2-component model-fit (§3) on these data, and find that the best-fit stellar and virial M/L ratios 
systematically vary from the values inferred from the original data (presented in §4) by ~ 2% and ~ 15% respectively, 
which are much smaller than the presented statistical errors. 

PSF Correction 

We next perform a PSF correction to remove any shape distortions of the background galaxies due to the tele- 
scope, leaving on l y thos e distortions caused by gravitational shear. Our correction scheme applies the method of 
iSchrabback et all ([2007D . which utilizes the KSB formalism along with modifications developed by iHoekstra et all 
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(fl99l and lErben etail IpOOl . 

While it should be possible to map the spatial variations of the PSF by directly measuring the shapes of the stars 
found in the data, we find that this is not the case for our lens sample, as there are too few stars per exposure to 
accurately sample the whole field. Instead, to increase our star count we use a series of well-sampled, low galactic 
latitude stellar field exposures to generate a series of polynomial PSF models. There is some concern in fitting a model 
PSF to our data instead of deriving it empirically, as the ACS PSF is temporally variable. However, this variability 
is primarily a function of a single pa r ameter: changes in the separation of the primary and secondary HST mirrors 
due to thermal breathing (|Kristll2003l iRhodes et aTll2007t ). By using a total of 181 exposures taken over the course of 



several months, our stellar-field PSF models span a wide range of mirror separations, allowing us to compile a nearly 
continuous database of varying PSF patterns that can be used to correct our science images. 

To create each model, we measure the anisotropy components (q*) and polarization tensor ratio (T*) - a parameter 
required for PSF correction that is related to the PSF width - as a function of position, using the stars in a given 
stellar field. The measured quantities are given by: 

5 * = (ps™*)-i e j ( A5 ) 
Tr [P sh *] 

where are the measured stellar ellipticity components, and P sh * and P sm * are the stars' shear and smear polarization 
tensors, respectively. 

Model stars are initially selected by simple cuts in size-magnitude space. A 0.8 pixel wide cut in half-light radius, 
centered on the stellar locus, is used to separate stars from galaxies. A bright-end magnitude cut is used to remove 
stars with saturated pixels and a SNR cut is used to remove stars that are too noisy to be accurately measured. From 
this initial star selection we generate 3rd-order polynomial models for the g ; * components, and a 5th-order polynomial 
model for T* . To improve the model, we evaluate each polynomial component at the coordinates of the stars, and 
compare the inferred values to the actual measured quantities. If, for any given model component, the observed value 
disagrees with the model by more than 2-ct, the star is considered an outlier (because of shape noise) and is removed 
from the sample. Once all outliers are removed, we regenerate the polynomial models with the remaining stars. This 
process is repeated iteratively until the number of stars remains stable, at which point we generate a final set of 
polynomials that represent the final PSF model for that stellar field. 

For each science exposure, we compare the measured PSF of each star in the field to the model PSFs, and determine 
the best-fit model by calculating the reduced \ 2 statistic given by: 

- [PSF*-PSF mod j (x i ,y i )] 2 



X) = E 1 ' : ^ (AT) 



where PSF* is the measured PSF of star i, and PSF mod j(x;, y;) is the model PSF of stellar template j, evaluated at 
the location (xj,yi) of the real star. 

Once best-fit PSF models are determined for all exposures associated with a composite science image, the coordinates 
of each model arc registered to the science image, and the models are combined by an exposure-time weighted average 
to create a final "master" PSF model used to correct each background galaxy, given by: 

E *i PSF i( x -y) A i 

PSF combo (x, y) = -i (A8) 

E ^ 

i 

where t; is the exposure time of a given exposure in the composite science image, and Aj is equal to 1 if the galaxy 
falls within the boundaries of the WFC chips in the exposure (not always true, due to dithering), and otherwise. 

As a test of our PSF correction, we plot the shapes of all stars found in our science images, both pre- and post- 
correction, (Figure [5]) ■ We do see an improvement: the stellar ellipticity components become more circular (ei )Pre = 
—0.0089, ei t p OS t = 0.0001; e2, pre = 0.0124, e 2iP ost = 0.0015), and show a mild reduction in spread (cr ei prs = 0.022, 
G &\ oat — 0-019; cr e2 r<! = 0.021,cr e2 oat = 0.016). By taking into account effects due to the PSF (and effects due to 
CTE degradation in the previous section) we are able to significantly reduce the dominant sources of systematic error 
involved in the weak lensing analysis. 

Data Cuts 

To improve the overall signal, we apply a series of data cuts to our master catalog. First, we exclude objects that 
have a half-light radius smaller than 3 pixels, as many of these objects are "bad" data objects such as stars or cosmic 
rays, which arc not affected by gravitational shear and only serve to dilute the overall shear signal. We also reject 
objects with a half-light radius larger than 14 pixels, and "bright" objects (itifsuw < 23), as these objects are likely 
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Figure 9. Result of PSF correction on all stars of our science images. Pre anisotropy-corrected shapes are plotted in red, while post- 
corrected shapes are shown in black. We can clearly see that the stars have become more circular after applying the combined PSF 
models. 

to be foreground interlopers. Additionally, we apply a conservative cut to galaxy shapes, rejecting any object with 
\gi\ > 2 (where the gi are the CTE + PSF corrected ellipticity components) as being unphysical or poorly shape- 
corrected in some way. Finally, we apply a general SNR cut to the data, rejecting any remaining object with an 
IMCAT significa nce (v) parameter le ss than 16; (as a comparison, this corresponds to an SNR < 4 cut using the SNR 
determination of lErben et al.l ()2001[ )'). After cutting the data, we are left with a total background galaxy density of 
~ 65 galaxies arcmin -2 (as compared to the ~ 100 galaxies arcmin -2 measured prior to cutting) which is high enough 
to measure the weak lensing signal around the galaxy stack. 



Catalog Stacking 

After correcting for PSF and CTE effects for obje c ts in an individual science field , we can apply t he standard weak 
lensing formalism (see e.g., iSeitz fc Schneider! 119971 ; iBartelmann fc Schneider! [2 00 It lWittmanTl2002j) to the field and 
determine an average mass profile from a stack of fields. 

To properly combine each lens field into a cohesive stack, we first apply a shift to the coordinates of each science 
field, such that the origin in each field corresponds to the location of the strong lens system. Because the lenses in 
our sample cover a range of redshifts, we convert each background galaxy's angular separation from the lens to a 
physical one, so that we can more easily stack galaxy fields together in catalog space. Finally, once all objects have 
been repositioned in terms of physical distances relative to a common center, we combine everything into a master 
galaxy catalog that contains information on galaxy position and corrected ellipticity, and information on the strong 
lens field from which it was initially drawn. In addition, geometric information about the strong lenses themselves are 
also included, which are critical in determining the physical mass profile at small scales. 
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Shear and Mass Profiles 

Since our master galaxy catalog contains information on both position and corrected cllipticity, we can combine 
these parameters to develop two radial shear profiles: a tangential "E-mode" profile given by: 



7t(R) = - 

and a 45° rotated "B-mode" profile, given by: 

7x (R) = -\ 



L gals,R 

E 



7i,i cos(20) + 72, i sin(20) 



-7i,i sin(20) + 72, j cos(20) 



(A9) 



(A10) 



where n ga is,R refers to the number of galaxies located within a radial bin centered at radius R and 71,, and 72,4 represent 
the shear components of galaxy i, equivalent to the components of reduced shear (c/i) in the limit of weak lcnsing. We 
calculate ~ft(R) and 7 X (R) in concentric circular annuli around the common center of our stacked lens sample. For a 
system with a single lens plane (or a stacked, coincident lens plane) the lensing signal measured with respect to the 
center of mass and averaged over all background galaxies should be contained entirely within the E-mode profile while 
the B-mode should vanish. Thus, by computing these two shear profiles, we can see information regarding the radial 
characteristics of the shear field from 7 t , as well as a check on systematic errors from 7 X . 

Next, we convert shear into a differential surface mass overdensity, AE, by scaling the tangential shear profile by 
the critical lensing density E cr it: 



^crit 



4ttG D,D la 



(All) 



a measurement of surface mass density that encodes the angular diameter distances between the observer and the lens 
(Di), the observer and the source (D s ), and the lens and the source (Di s ). The AE(ii) profile is thus given by: 



AS(i?) = E(< R) - E(R) = £ cri t7t(#) 



(A12) 



where E(< R) represents the average surface mass density enclosed within projected radius R. Since E cr ;t is a function 
of both lens and source redshifts (zi and z s , respectively), in principle we must construct a unique T,„i t {zi, z s ) for every 
source galaxy that is used in the mass overdensity profile. In practice, however, this is impractical since a large number 
of background galaxies do not have measured redshifts. Instead, as in G07, we define an average critical density for 
each lens field: 



\4irGJ 



1 



(A13) 



where w(zi) is a weighted factor of (D s /Di s 
is given by: 



w(zi) 



Diw(z{] 

designed to fully account for background source redshift distribution. It 

/ dn(z s )\ D is (z s ) 



dz s 



V dz s 



D s (z s ) 



(A14) 



where dn(z s ) / dz s is the redshift distribution of background source galaxies. Since our image exposure times and data 
cuts are similar to G07 we use their redshift distribution, extrapolated from COSMOS: 



dn(z s ) 
dz„ 



jr(o/6) 



-^/^\ Zs /z r 



(A15) 



with zq = 0.345, a = 3.89, and 6 = 1. From these parameters, we calculate a mean background source redshift of 
= 1.34. A full list of w(zi) parameters can be seen in Table[TJ 
We next compute an estimator of AE(i?): 



AE(iJ) = 





JV sourco ,, 
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E crit,J E 7 *>' CT e 
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(A16) 



where iV sourcc .; is the number of background galaxies associated with lens field I, E^ rit ; is the average critical density 
for lens field /, and o~ e j is the uncert ainty assigned to t h e tan gential shear estimate of galaxy i, which we take to 
be the shape uncertainty as defined in lGavazzi fc Soucaill (|2007f ) to be proportional to a galaxy's shape measurement 
error added in quadrature with the intrinsic shape noise of all galaxies (ctq = 0.3) to prevent over-weighting. With this 
conversion, we are able to transform shear to mass, and develop a mass overdensity profile for our data. 
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As a check on the systematics of our assumed redshift distribution, we vary the parameters of the distribution, 
calculating a unique set of {a, b, zo} parameters for each lens field. We do this by measuring the observed F814W 
magnitude distribution of the galax ies in the field and convert this to a redshift distribution, using the method 
described in iSchrabback et al.l (|2007f ). After doing this, we recalculate the SJ. rit , mJ(zz), and AS(i?) estimator values 
for the entire sample, and repeat the weak lensing analysis with these updated values. We find that, regardless of the 
choice of redshift distribution parameters, the best-fit model parameters determined here (as well as those described 
in §3) vary by < 7%. Since this is smaller than the statistical uncertainty of our measurements, we conclude that 
uncertainties in the redshift distribution are negligible for our purposes. 



